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Confirmation of the Copernican principle at Gpc radial scale and above from the 
kinetic Sunyaev Zel'dovich effect power spectrum 
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The Copernican principle, a cornerstone of modern cosmology, remains largely unproven at Gpc 
radial scale and above. Here we will show that, violations of this type will inevitably cause a first 
order anisotropic kinetic Sunyaev Zel'dovich (kSZ) effect. If large scale radial inhomogeneities have 
amplitude large enough to explain the "dark energy" phenomena, the induced kSZ power spectrum 
will be much larger than the ACT/SPT upper limit. This single test confirms the Copernican 
principle and rules out the adiabatic void model as a viable alternative to dark energy. 
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Introduction. — The Copernican principle has been 
a fundamental tenet of modern science since the 16" 1 
century and is also a cornerstone of modern cosmology. 
It states that we should not live in a special region of 
the universe. Cosmic microwave background (CMB) ob- 
servations verify the statistical homogeneity of the last 
scattering surface fl. Galaxy surveys verify the radial 
homogeneity up to the Gpc scale [2j. However, radial 
homogeneity at larger scales remains unproven. 

Testing the Copernican principle is of crucial impor- 
tance for fundamental cosmology. If the Copernican prin- 
ciple is violated such that we live in or near the center of 
a large (~ Gpc) void as described by a Lematre-Tolman- 
Bondi (LTB) space-time |3| in which the matter distribu- 
tion is spherically symmetric, the apparent cosmic accel- 
eration i], [B[ can be explained without cosmological con- 
stant, dark energy or modifications of general relativity 
Throughout the paper we will restrict to this type 
of violation of the Copernican principle. Various tests 
of the Copernican principle have been proposed and a 
large class of void models has been ruled out (e.g. 043)- 
Here we propose a powerful single test which confirms 
the Copernican principle at Gpc radial scale. 

The kSZ test. — A generic consequence of violating 
the Copernican principle is that some regions will ex- 
pand faster or slower than others and as photons tran- 
sit between these regions there will be a relative motion 
between the average matter frame and CMB. When rel- 
ative motions between free electrons and photons exist 
the inverse Compton scattering will induce a shift of the 
brightness temperature of CMB photons via the kinetic 
Sunyaev Zel'dovich (kSZ) effect [10(. This temperature 
shift will be anisotropic on our sky tracing the anisotropy 
of the projected free electron surface density. This test 
of the Copernican principle has been applied to cluster 
kSZ observations [3,[§|, where the electron surface density 
is high. However this effect applies to all free electrons 
which exist in great abundance everywhere in the uni- 
verse up to the reionization epoch at redshift z > 6 (and 
comoving distance ^ 6/i _1 Gpc), whereas clusters are rare 



above z ~ 1. So one can expect a more sensitive test 
from blank field CMB anisotropy power spectrum mea- 
surements than from cluster measurements as has been 
demonstrated for the "dark flow" |ll| induced small scale 
kSZ effect 0. 

Free electrons have local motion vl with respect to 
the average matter frame and the subscript "L" refers 
to "local" . It vanishes when averaging over sufficiently 
large scale. However, when the Copernican principle is 
violated at large scale, electrons will have relative mo- 
tion vn between the average matter frame and CMB, 
which does not vanish even when averaging over the Hub- 
ble scale. Correspondingly the induced kSZ temperature 
fluctuation T3, 12\ has two contributions, 



AT(h) = AT L {h) + AT H (n) 



(1) 



The first term on the r.h.s is the conventional kSZ effect, 
AT L (h) = Tcmb x f[l + 5 e (n, z)] ^ ( "' z) ' H dr e . (2) 



Here, n is the radial direction on the sky. r c is the mean 
Thomson optical depth to the corresponding redshift and 
S e is the fractional fluctuation in the free electron number 
density. The last term in Eq. [I] is new and does not 
vanish in a non-Copernican universe, 



AT H (n) = T CMB x l[l+6 e (n,z)} VR( - n ' z) ' n dT e 



9.1/xK 



v-g ■ n 8 c {n, z) dr e 
10 4 km/s 0.1 0.001 



(3) 



The last expression neglects the J wh • ndr e term, which 
has no direction dependence in LTB models in which 
we live at the center, and is therefore not observable. 

varies slowly along radial direction and does not suf- 
fer the cancellation of vt, in the conventional kSZ effect 
T3I 14 1 . The small scale anisotropy power spectrum will 
be quadratic in the amplitude of S e (which does fluctuate 
about zero) so we can say that ATh/T is first order in 
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the density fluctuations. Throughout this paper, unless 
otherwise specified, we will focus on this linear kSZ ef- 
fect. We restrict ourselves to adiabatic voids in which 
the initial matter, radiation, and baryon densities track 
each other. This is what one would expect if baryoge- 
nesis and dark matter decoupling occurs after the pro- 
cess which generates the void inhomogeneity. We also 
restrict ourselves to voids outside of which both matter 
and radiation are homogeneous. Adding additional inho- 
mogeneities will generically lead to larger values of vh- 

To explain the dimming of SNe-Ia and hence the ap- 
parent cosmic acceleration without dark energy and mod- 
ifications of general relativity, we shall live in an under- 
dense region (void) of size > l/i _1 Gpc, with a typical 
outward velocity «h ^ 10 4 km/s (e.g. |9|). Given the 
baryon density ft^h 2 = 0.02 ± 0.002 from the big bang 
nucleosynthesis [l5|, r e > 10~ 3 . Scaling the observed 
weak lensing rms convergence K^10~ 2 at^7 17], the 
rms fluctuation in 6 C projected over Gpc length is > 0.1 
at such scale. Hence such a void generates a kSZ power 
spectrum AT- 2 - > 80/xk 2 at multipole I = 3000. This is 
in conflict with recent kSZ observations. The South Pole 
telescope (SPT) collaboration [HI found AT 2 < 6.5/xK 2 
(95% upper limit) and the Atacama cosmology telescope 
(ACT) collaboration [3| found AT 2 < 8^K 2 . This sim- 
ple order of magnitude estimation demonstrates the po- 
tential discriminating power of the kSZ power spectrum 
measurement. It suggests that a wide range of void mod- 
els capable of replacing dark energy are ruled out. This 
also demonstrates how purely empirical measurements 
of CMB anisotropics and the large scale structure (e.g. 
weak lensing) can in principle be combined to limit non- 
Copernican models without any assumptions of how the 
inhomogeneities vary with distance. 

We perform quantitative calculation for a popular void 
model, namely the Hubble bubble model (|8| and refer- 
ences therein). In this model, we live at the center of a 
Hubble bubble of constant matter density J7o < 1 embed- 
ded in a flat Einstein-de Sitter universe (fi m = 1). The 
void extends to redshift z dg , surrounded by a compen- 
sating shell (zcdgo < z < z ut) and then the flat Einstein- 
de Sitter universe (z > z ou t)- The kSZ effect in this uni- 
verse has two components, (1) the linear kSZ arising from 
the large angular scale anisotropies generated by matter 
(a) inside the void, (b) in the compensating shell, (c) out- 
side the void; (2) the conventional kSZ effect quadratic 
in density fluctuation [l4| and the kSZ effect from patchy 
reionization [2]]]. The contributions of each of these to 
the anisotropy power spectrum are uncorrelated. Hence 
the ACT/SPT measurements put an up per limit on the 
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total. The later contributes ~ 3.5//K 2 
left for the first component is < 3/jK 2 . However, we 
will test the Copernican principle in a conservative way, 
by requiring the power spectrum of the first component 
generated by matter inside the void to be below the SPT 
upper limit 6.5/zK 2 at I = 3000. 



For a general Hubble bubble «h is determined by both 
Doppler and Sachs- Wolfe anisotropies generated by the 
void and depends qualitatively on the size of the void 
[H [23[ . As we shall see below it is only Hubble bubbles 
with z e( jgc < 1 which are consistent with both the SNe 
data and the proposed kSZ test, and for these a simple 
Doppler formula can be used @, Hi] 



vh(z) « [Hi{z)-H e ] 



A&,co(z) 



1 



(4) 



where, Hi(z) is the Hubble expansion rate inside the void 
as a function of redshift, H c gives the Hubble expansion 
rate exterior to the void at the same cosmological time, 
D\ iCO {z) is the comoving angular diameter distance to 
redshift z. 

The temperature fluctuation at multipole t generated 
by the linear kSZ effect inside of the Hubble bubble is, 
under the Limber approximation, 



AT 2 {I) = (9.1 M K) 2 
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Da,co(z) 
c/Hi(z) 



(5) 



dz . 



Here AT 2 = T 2 MB £(£ 



l)Ci/(2n), Ci is the correspond- 
ing angular power spectrum, A 2 (fc,z) = £^P e (k,z) 
is the dimensionless electron number overdensity at 
wavenumber fc and redshift z. 

In our calculations we approximate P e by the matter 
power spectrum P m and approximate P m by it's form in a 
standard ACDM cosmology. It is non-trivial to calculate 
P m in LTB models, since even at linear scales the ex- 
pansion rate is locally anisotropic so the inhomogeneities 
will have an anisotropic power spectrum (see 16J): and 
since we are no longer assuming the cosmological prin- 
ciple one could also expect large scale variations in the 
initial inhomogeneities. The measured matter clustering 
and its evolution agree with the standard ACDM cosmol- 
ogy to a factor of ~ 2 uncertainty up to z ~ 1 [17J, [25j , 



as do the galaxy clustering and evolution 26|. A mini- 
malist approach is to simply use the ACDM predictions 
since any viable LTB models must be consistent with 
this data. If this assumption are not satisfied then one 
should be able to obtain even tighter constraint by con- 
sidering these extra tests. Here we use P m calculated by 
the CMBFAST package [27j j . nonlinear clustering from 
the halofit formula [28[ all using assuming ACDM with 
n m = 0.27, n A = 1 - fi m , n h = 0.044, a s = 0.84 and 
h = 0.71. All other quantities such as r c and vjj are cal- 
culated based on the void model with the same fib and 
Hi(z = 0) = 100ft, km/s/Mpc. The kSZ power spectrum 
is then computed using Eq. [5l 

Constraints on the void model. The ACT/SPT 
upper limit rules out large voids with low density (Fig. 
[T]). Only those voids either with f2 -t 1 (£Iq > 0.8) 
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FIG. 1: The kSZ test. Black curves have constant ATr(£ = 
3000). The thick one highlights the SPT 95% upper limit, 
AT 2 < 6.5/iK 2 PH. The kSZ test alone rules out large voids 
with low density and strongly supports the Copernican prin- 
ciple. The dashed and dotted contours are the 2-cr and 3-ct 
constraints from the UNION2 supernova data [3(|. The kSZ 
test robustly excludes the Hubble bubble model as a viable 
alternative to dark energy. 



or Zcdgo — > (zodgc ^) 0.2, corresponding to void radius 
< 0.6/i _1 Gpc) survive this test (Fig. [1]). These results 
agree fairly well with those of in a more recent paper (Fig. 
6, [111), who have used a more sophisticated treatment. 1 
The kSZ test is highly complementary to other tests 
such as the supernova test. Our SNe la constraint follows 
rcf. 8] but uses the improved UNION2 data with 557 SNe 
la [30|. Not allowing for additional intrinsic dispersion 
of the SNe magnitudes we find a minimum x 2 is 605. 4. 2 
Hubble bubble models within 3er contour have typical 
AT 2 > 10 3 ^K 2 at I = 3000, two orders of magnitude 
larger than the SPT upper limit 6.5/xK 2 [l8j]- On the 
other hand, Hubble bubble models consistent with the 
SPT result have A X 2 > 209 {\ 2 > 814) for the SN la 
test and hence fail too. Thus the combination of SN la 
observations with small scale CMB anisotropy apparently 
rule out all Hubble bubble models. 

Our kSZ calculation is based on these assumptions: 1) 



1 However, since [29j uses a different smooth void model, our re- 
sults are not directly comparable. 

2 Although this indicates a poor fit including systematic errors and 
intrinsic magnitude dispersions would improve the fit. 



FIG. 2: The maximal deviation from the overall expansion 
allowed by the SPT observation, for each mass shell of Az = 
0.4, which corresponds to l/i _1 Gpc at z ~ 0, 0.7/i -1 Gpc at 
z ~ 1 and 0.5h~ Gpc at ~ 2. 



Ob/i 2 is the same as in the standard BBN analysis, 2) P m 
is the same as in a ACDM model, (based on the argument 
that any viable void model must reproduce the observed 
matter clustering), 3) P c = P m (good to ~ 10% accuracy 
f20(), 4) eq. [j] fo r velocities (roughly accurate for sub- 
horizon voids [24j which is required by CMB data [1,0]), 
5) neglect of kSZ contributions from the compensating 
shell (which would only increase kSZ anisotropy), 6) a 
simple adiabatic Hubble bubble void, and 7) no CMB 
flow (intrinsic dipole) from non-adiabatic initial condi- 
tions outside the void. We expect that relaxing l)-6) in 
reasonable ways could not significantly reduce the ten- 
sion imposed by the kSZ test, since for void models to 
explain the observed SN dimming, they must have large 
scale gravitational potential of large amplitude and hence 
must have large vh and large kSZ effect. For example, 
{29j adopted a void model of different profile and found 
much weaker SN constraint, but the generated kSZ power 
is nevertheless much larger than the ACT/SPT upper 
limit. This demonstrates the great discriminating power 
of the kSZ test. Completely relaxing 7) could change our 
conclusion for rather contrived initial conditions [31| , but 
would generically lead to even larger and more unaccept- 
able kSZ effect. Thus comparing kSZ with SNe is by 
far the most stringent test of the void models and the 
Copernican principle at Gpc scale and above. We con- 
clude that any adiabatic void models capable of explain- 
ing the supernova Hubble diagram would likely generate 
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too much kSZ power on the sky to be consistent with the 
ACT/SPT upper limit. This strengthens the evidences 
for cosmic acceleration and dark energy. 

Constraints on the Hubble flow. — Still, violation 
of the Copernican princi ple less dramatic than the above 
void models may exist [32j [. For example, there could 
be large scale density modulation on the ACDM back- 
ground. As long as the amplitude of the modulation is 
sufficiently small, it can pass the supernova test and the 
structure growth rate test. However, if unaccounted, it 
could bias the dark energy constraint. The kSZ test is 
able to put interesting constraint on this type of viola- 
tion. We take a model independent approach and pa- 
rameterize the violation of the Copernican principle by 
AH(z), the deviation of the Hubble expansion of a mass 
shell of size Az centered at redshift z from the overall 
expansion of the background universe. The ACT result 
constrains \AH(z)/H(z)\ < 1% for each mass shell of ra- 
dial width ~ l/i _1 Gpc (Fig. [5]). This estimation neglect 
contributions from other mass shells so the actual con- 
straint is tighter. This test can be carried out on each 
patch of the sky to test the isotropy of the Hubble flow. 

The above test is not able to determine at which red- 
shift a violation of the Copernican principle occurs, since 
the kSZ power spectrum is the sum over all contribu- 
tions along the line-of-sight and hence has no redshift 
information. This problem can be solved with the aid of 
a survey of the large scale structure (LSS) with redshift 
information. 

The basic idea is the same as the one proposed by 
[T3 | to probe the dark flow through the kSZ-LSS density 
distribution two point cross correlation. This cross cor- 
relation is non-zero only in non- Copernican Universes, 
since the velocity vh varies slowly over the clustering 
length of the LSS and since the linear kSZ effect is lin- 
ear in density. Since the cross correlation vanishes for 
the conventional kSZ effect, a non- vanishing cross corre- 
lation signal can serve as a smoking gun of violation of 
the Copernican principle. The thermal SZ contaminates 
the measurement. However, it can be largely removed by 
spectral fitting or observing at its null: 217 GHz. Since 
the redshift surveys can map the LSS with much higher 
S/N than kSZ measurements, this cross correlation can 
achieve much higher S/N than the kSZ auto-correlations. 
We thus expect that small scale CMB anisotropy surveys, 
such as ACT and SPT, in combination with deep LSS 
surveys will be able to put more stringent constraints on 
violations of the Copernican principle at each redshift 
and each direction of the sky. 
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